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Abstract. We investigate the sensitivity of the Gamma-ray Large Area Space 
Telescope (GLAST) to indirectly detect weakly interacting massive particles (WIMPs) 
through the 7-ray signal that their pair annihilation produces. WIMPs are among the 
favorite candidates to explain the compelling evidence that about 80% of the mass in 
the Universe is non-baryonic dark matter (DM). They are serendipitously motivated by 
various extensions of the standard model of particle physics such as Supersymmetry and 
Universal Extra Dimensions (UED). With its unprecedented sensitivity and its very 
large energy range (20 MeV to more than 300 GeV) the main instrument on board the 
GLAST satellite, the Large Area Telescope (LAT) , will open a new window of discovery. 
As our estimates show, the LAT will be able to detect an indirect DM signature for 
a large class of WIMP models given a cuspy profile for the DM distribution. Using 
the current state of the art Monte Carlo and event reconstruction software developed 
within the LAT collaboration, we present preliminary sensitivity studies for several 
possible sources inside and outside the Galaxy. We also discuss the potential of the 
LAT to detect UED via the electron/positron channel. Diffuse background modeling 
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and other background issues that will be important in setting limits or seeing a signal 
are presented. 
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1. Introduction 

The Gamma-ray Large Area Space Telescope (GLAST) [H El [3l H] is a satellite-borne 
7-ray detector launched on 11 June 2008. The Large Area Telescope (LAT) is the main 
instrument on GLAST, which will also host the Gamma-ray Burst Monitor (GBM) 0. 
In this paper, we investigate the potential of the LAT for confirming, or constraining, 
the most interesting models of particle dark matter of the Universe, WIMP (weakly 
interacting massive particle) models. The model of the lightest supersymmetric particle 
is the most studied template, which we use for most of our analysis. For general reviews 
on supersymmetric and other models of particle candidates for dark matter, see [6]. We 
also study some aspects of a completely different class of models, so-called Kaluza-Klein 
models of Universal Extra Dimensions (see [7] and references therein). Although esti- 
mates of 7-ray signals exist in the literature (see [6] ) , the new feature in the present paper 
is that the theoretical predictions are fed through the experimental response function of 
the LAT and the software analysis chain that will be used for the actual data analysis 
after launch. Thus we give good estimates for the potential of GLAST for detecting 
or limiting dark matter models. We also discuss the most important astrophysical and 
instrumental backgrounds. 

The LAT will have more than an order of magnitude better sensitivity in the 20 MeV 
to 10 GeV region than its predecessor, EGRET onboard the Compton Gamma-ray Ob- 
servatory [8], and furthermore will extend the high energy region to roughly 300 GeV. 
Therefore, the LAT emerges as an instrument that is well suited to search for signals 
from dark matter annihilation, which in the case of WIMPs should be populate just this 
energy range. As the Large Hadron Collider (LHC) at CERN will also start taking data 
by the end of 2008, there is a non negligible probability to detect a good WIMP candi- 
date, and verify through the 7-ray signal that such particles constitute the dark matter 
halo of the Milky Way (or neighboring galaxies or sub-halos, which we also treat). For 
a thorough discussion of the interplay between discovery at LHC and detection of dark 
matter through other methods, in particular 7-ray observations, see [9]. 

It is now established beyond reasonable doubt from a combined study of the cosmic mi- 
crowave background radiation |10j, supernova cosmology [llj and large Galaxy redshift 
surveys fl2[ [13] that non-baryonic dark matter is needed, while other models such as 
modifications of the laws of gravity have problems explaining the wealth of observations 
that now are in place. For example, it seems that the combined X-ray and optical ob- 
servations of the "Bullet Cluster" essentially exclude explanations not involving dark 
matter [Ti] . 

I For a summary of overall design properties, general performance, planned operation modes, data 
analysis tools and strategies of the LAT, as well as a full list of participating scientists, see the web site 
of the collaboration [5]. 



4 



Technically denoted CDM, Cold Dark Matter, the particles constituting the cosmo- 
logically required dark matter have to be moving non-relativistically at the epoch of 
structure formation to reproduce the observed structure of the Universe, especially at 
small scales. This property is always fulfilled by particles with masses in the GeV to 
TeV range that interact with the weak interaction strength, i.e. WIMPs. They will 
have velocities which have redshifted since the time of thermal decoupling in the early 
Universe, and they will now move with typical Galactic velocities f/c ~ 10^'^ in the 
Milky Way halo. This is in contrast with massive neutrinos. Hot Dark Matter, which 
give an observationally disfavored top-down structure formation scenario with relatively 
little structure on small scales. In fact, the agreement of the cosmological power spec- 
trum with that of CDM allows us to put stringent bounds on the neutrino mass (see, 
e.g, [IS]). On the other hand, the non-zero neutrino masses inferred from measured 
neutrino oscillations, although not enough to explain more than a few percent of the 
dark matter, constitute a first demonstration that non-baryonic dark matter indeed is 
likely to exist. 

We thus have an excellent class of particle candidates, WIMPs, which behaves as the 
cosmologically needed CDM, and which could have been thermally produced in the early 
Universe to give the required relic density [10] VLcDuh'^ ~ 0.1, where h is the Hubble 
constant in units of 100 kms~^Mpc~^. The ability of WIMPs to give the measured relic 
density from readily computed thermal processes in the early Universe without much 
fine tuning is sometimes termed the 'WIMP miracle" . 

Of course it has to be kept in mind that dark matter does not necessarily have to be 
WIMPs in the mass range detectable by GLAST. The Warm Dark Matter model with 
thermal-relic particle masses above 2 keV may also explain the present observations in 
both sky surveys and N-body simulations. In the case when, say, a gravitino is the 
LSP, no signal would be observable with GLAST. Here we work, however, with the 
assumption that the dark matter particle has detectable 7-ray couplings and present 
the discovery potential for GLAST. 

In Figure 1, we sketch how 7-rays are produced from DM annihilation, and also show 
a schematic of a typical simulation and analysis chain as used in this paper. Table [1] 
shows the various approaches to the search for WIMP dark matter signals in 7-ray data 
explored in this paper. A "smoking gun" would be the detection of line emission in 
WIMP annihilation [16], through the loop-induced annihilation into two photons, which 
for slow-moving dark matter particles would give rise to a striking, almost monoen- 
ergetic, photon signal. However, the branching ratio for the annihilation into lines is 
typically about 10~^ or less in most models, as WIMPs turn out to be electrically neu- 
tral, and thus do not couple directly to photons. There are, however, some exceptions to 
this estimate [IT]. In addition to considering signatures in 7-rays we also illustrate the 
LAT capability to detect electron/positrons, which could provide signatures of Kaluza- 
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Table 1. The various venues GLAST will explore in its search for WIMPs, and the 
advantage/disadvantage of each method. 
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Klein particle Dark Matter. 

The paper is organized as follows: In section [2] we give a description of the LAT instru- 
ment and the software used for the analyses presented in this paper. In section [3] the 
calculations of the WIMP signal flux are described and a discussion of the considered 
7-ray background is given. In section H] we summarize the sensitivities to generic WIMP 
annihilation signals achievable by the LAT for the search channels currently pursued by 
the LAT collaboration. Specific models are studied in section |5l 

2. The LAT 

The main instrument on board GLAST is the LAT, which is complemented by the 
GBM. The LAT is modular, consisting of a 4 x 4 array of identical towers. Each 40 x 
40 cm^ tower is composed of a tracker, a calorimeter, and a data acquisition module. 
The tracker array is covered by a segmented anti-coincidence shield (ACD) . The tracking 
detector consists of 18 xy layers of silicon strip detectors interleaved with 16 layers of 
tungsten foils. 

Each calorimeter module has 96 CsI(Tl) crystals, arranged in an 8 layer hodoscopic 
configuration with a total depth of 8.6 radiation lengths, giving both longitudinal and 
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Figure 1. A diagrammatic flow of how gamma rays are produced by annihilation 
of dark matter and elements of the analysis chain used by the GLAST collaboration 
to detect them. The double question mark in the simulation chain indicates high 
uncertainty in the models of dark matter density and the new particle theories 
discussed in the paper. The single question mark over the cosmic ray propagation and 
interaction models indicates lesser, although significant, uncertainty in those models 
that generate backgrounds to the potential dark matter gamma ray signal. In this 
paper GALPROP ('section r3.2l) is used to estimate those backgrounds. In the next step, 
7-ray detection is simulated using standard detector simulation packages (GEANT 4). 
Finally,these simulated LAT events are treated by various analysis software programs 
(event reconstruction and statistical analysis) to generate the results presented in this 
work. The same procedure is applied to the smoking gun signal of XX ~^ 11^ except 
that in this case hadronlzation does not have to be taken into account. 



transverse information about the energy deposition pattern The calorimeter's depth 
and segmentation enable the high-energy reach of the LAT and contribute significantly 
to background rejection. The ACD is the LAT's first line of defense against the charged 
cosmic ray background. It consists of 89 different size plastic scintillator tiles and 
9 ribbons with wave-length shifting fiber readout. The segmentation is necessary to 
suppress self-veto effects caused by secondary particles emanating from the calorimeter 
showers of high energy 7-rays [18] . 

2. 1 . LA T Exposure 

For this paper, simulations of LAT all-sky "exposures" of 2 months, 1 year, 5 years 
and 10 years are used in the analyses. LAT exposure is defined as the amount of cm^ 
s the LAT effective area integrates over many orbits, which is a complex calculation. 

§ With the tracker the LAT presents 10 radiation lengths for normal incidence. 



2.0E+11 2Jf:-Ht1 2.3C-!-11 Z.4E+-11 

Figure 2. The simulated LAT exposure for 5 years of all-sky scan. The effect of 
turning off the LAT while in the SAA is included in the exposure. This exposure is 
calculated for a photon energy of 100 GeV. The plot is in Galactic coordinates with 
the values of exposure shown on the grey (coloured in coloured versions of the paper) 
bar in units of cm^ sec. 



The GLAST orbit will be oriented at 25.4 deg. to the Earth's equator at an altitude 
of 565 km. Except where otherwise noted, all the analyses in this paper use codes that 
calculate the exposure from all-sky scanning consisting of alternating orbits with the 
LAT rocked by 35 deg toward the north celestial pole and then the south pole relative to 
local zenith. The GLAST orbit has about a 52.5-day precession period over which this 
all-sky scan achieves better than the required ± 20% uniformity of exposure on the sky 
(not uniformity of sensitivity due to backgrounds that vary over the sky). Thus all of the 
analyses presented in this paper have the required uniformity of exposure since they all 
integrate longer than 52.5 days, but it should be noted that the only analyses making use 
of this fact are the line search (section 14.51) . the search for Galactic satellites (section 
14. 3p and the sensitivity to mSUGRA models (section [5]). The other analyses correct 
for the effects of exposure. Figure [2] shows the simulated LAT exposure in Galactic 
coordinates, for 5 years of all-sky scan, in units of cm^s. The effect of turning off the 
LAT while in the South Atlantic Anomaly (SAA) is included in the exposure. This 
exposure is calculated for a photon energy of 100 GeV and is insensitive to the energy 
over the energy range of interest in this paper. The simulated exposure is uniform over 
the sky to ± 15 %. 



8 



2.2. LAT Analysis Software 

The simulations and analyses presented in this paper make use of the LAT Monte Carlo 
simulation, GLEAM, and the Science Tools software package. GLEAM is a C++ package 
based on the GEANT 4 toolkit. GLEAM provides a detailed simulation of the physics 
processes and digitization within the LAT using a highly detailed geometrical description 
of the detector. It has been validated in extensive beam tests. The ScienceTools [19] 
are being developed jointly by the LAT collaboration and the GLAST Science Support 
Center. The ScienceTools package comprises astronomical analysis software at a level 
similar to that available for other high energy satellite missions, such as Chandra or 
XMM, and these tools will be distributed in a similar manner for guest observers to 
perform analyses of LAT data. The two main programs used in the current paper 
are the observation simulation tool gtobssim and the source analysis tool gtlike. 
The observation simulation tool allows users to define celestial 7-ray sources with any 
spectral properties. These sources can be point-like or spatially extended. For the 
latter, the angular distribution of the emission is given using FITS image files. Given 
incident photons generated according to the input source models, gtobssim uses the 
LAT instrument response functions (IRFs) to select photons as detected in the LAT 
and to assign measured energies and directions. The IRFs, i.e., the effective area, energy 
dispersion, and point spread function, are ascertained from the detailed instrument 
Monte Carlo simulations. The source analysis tool, gtlike, uses a maximum likelihood 
method to fit source spectral parameters such as flux and power-law spectral index, 
though more complex spectral models are available. Since the detected counts for sources 
near the detection limit will be fairly low, gtlike calculates a likelihood function based 
on the Poisson probability using the source model folded through the LAT IRFs [20| to 
provide the expected model count. 

The process of the optimization of the event selection for high level analysis is 
ongoing. In the analyses presented in this paper, the event data are binned in direction 
and energy and gtlike optimizes the binned log-likelihood. A number of ancillary 
ScienceTools are also used that create the counts and exposure maps for the gtlike 
program. 

3. Estimating the expected signal and background 

In this section we give a general description of how we estimate the sensitivity of the 
LAT to 7-rays originating from WIMP annihilation. We describe the general way of 
estimating the signal flux and the backgrounds we consider in most of our searches. 
It should be noted, that in all our analyses we assume that the WIMP annihilation 
spectrum is the same no matter where observed (except for the redshift effect for the 
cosmological WIMP annihilation case), i.e. we assume the astrophysical environment 
does not affect the spectrum. 
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After calculating WIMP annihilation and background fluxes, we use the fast GLAST 
detector simulation as described in the previous section for simulating LAT data. For 
the estimate of the sensitivity to specific particle physics models (see section [5]), we 
assume a constant exposure to simplify computations. Finally, the sensitivity is calcu- 
lated employing variants of the above mentioned maximum likelihood method unless 
otherwise stated. 



To begin with, we study the LAT sensitivity to generic WIMP annihilation models, 
i.e. we simply consider the annihilation of a WIMP without any specific predictions 
on branching ratios and annihilation cross section, which could be provided by some 
underlying theory. This choice is motivated by lack of agreement between numerical 
computations. Various studies (see [2ll |22] and references therein) show, that in the 
framework of the widely popular and most restrictive supersymmetric extension of the 
Standard Model, the minimal supergravity framework (mSUGRA), there are important 
differences in the outputs of different spectrum calculators and Renormalization Group 
Evolution codes. This is particularly true at large values of tan{(3), the ratio of vacuum 
expectation values of the two neutral component of the SU (2) Higgs doublet (see sec- 
tion [5]), or for large values of mo, the common scalar masses. The practical consequence 
of these differences is uncertainty in the predictions of the various experimental observ- 
ables like relic density or annihilation cross-sections. 



Nevertheless, especially in the context of comparison with accelerator experiments, it 
is interesting to consider specific particle physics models. Two examples are studied in 
section [5] 



3.1. Calculation of the signal -ray flux 

The 7-ray continuum flux from WIMP annihilation at a given photon energy E from a 
direction that forms an angle between the direction of the Galactic center and that 
of observation is given by [23] 



The particle physics model enters through the WIMP mass my/ imp-, the total mean an- 
nihilation cross-section a multiplied by the relative velocity of the particles (in the limit 
of f — > 0), and the sum of all the photon yields dNf/dE for each annihilation channel 
weighted by the corresponding branching ratio Bf. As pointed out in [21], apart from 
the T~^T~ channel, the photon yields are quite similar. 

The integral in Equation ([1]) is the integral along the line of sight (l.o.s) of the assumed 
density squared, p(/)^, of WIMPs. Among the kinematically allowed tree-level final 
states, the leading channels are often bb,tt,T~^T~ ,W^W~ , Z^Z^. This is the case, 
e.g., for neutralinos and, more generically, for any Majorana fermion WIMP. For such 
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particles the S-wave annihilation rate into the light fermion species is suppressed by 
the factor rn'j/m'^jj^p, where is the mass of the fermion in the final state. The 
fragmentation and/or the decay of the tree-level annihilation states gives rise to photons. 
The dominant intermediate step is the generation of neutral pions and their decay into 
27. The simulation of the photon yield is standard. We take advantage of a simulation 
performed with the Lund Monte Carlo program Pythia 6.202 [25] implemented in the 
DarkSUSY package p3]. The density distribution in DM halos, from simulations, are 
well fitted by simple analytical forms, where the most common one is given by 

P(^) = T^^- (2) 



This function behaves approximately as a broken power-law that scales as r~'^ close 
to the center of the halo, r~^^~'^^^°' at an intermediate distance r^, and in the 
outskirts of the halo. In the remainder of the paper, we usually assume (a, /3, 7,rs) = 
(1, 3, 1,20 kpc), corresponding to the Navarro-Frenk- White (NFW) [26] profile, and 
(a, P, 7,rs) = (1, 3, 1.5,28 kpc), corresponding to the Moore profile [27]. Unless 
otherwise stated, we normalize ps to give 0.3 GeV cm~^ at the Sun and our distance to 
the Galactic center is assumed to be 8 kpc p8] . These kinds of profiles increase sharply 
towards the origin and are therefore called 'cuspy'. Non-divergent profiles, referred to 
as cored, have also been proposed. We also give estimates for cored profiles, where 
applicable. 

For the case of the cosmological WIMP annihilation, the yield of 7-rays has to be 
redshifted and the line of sight integral has to be performed over halos of all redshifts, 
which implies modeling their structure as a function of redshift. For more details see 
section 14. 6[ 



3.2. Backgrounds 

The sensitivities to a DM signal presented in this paper depend critically on accurate 
estimates of the following backgrounds: diffuse Galactic 7-rays, extragalactic diffuse 
7-rays, and charged particles in the instrument. In this section we discuss these dif- 
ferent contributions and how we use them in the estimate of the sensitivity of GLAST 
to a DM signal. As we will see, the prediction of the conventional 7-ray diffuse emis- 
sion is subject to considerable freedom. A quantitative treatment of the systematic 
uncertainties introduced in this prediction requires mapping a likelihood function in a 
multi-dimensional parameter space with non-linear dependences. We assume that the 
diffuse emission models that we use is the true background and the evaluation of the 
background systematics is beyond the scope of this paper. 



We also do not take into account the contribution of unresolved point sources. The 
contribution of unresolved point sources is small compared with the true Galactic in- 
terstellar emission. The number and flux distribution of EGRET point sources implies 
that the contribution of unresolved sources to the Galactic diffuse emission is about 10 
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% or less. The contribution will be evensmaller in a LAT measurement because fainter 
sources will be resolved. The contribution of point sources above the EGRET energy 
range is poorly constrained due to lack of observations, but is also likely to be small [29] 

For this reason (and considering the significant astrophysical uncertainties which enter 
the calculation of signal fluxes) we content ourselves with two illustrative examples, 
chosen to give a rough idea of what the systematic uncertainties might be: the 
"conventional" GAL PROP [30J and the "optimized" GAL PROP model [31] , which are 
discussed in detail below. For more details on the GALPROP code, see [32]. For recent 
reviews on cosmic ray (CR) propagation and diffuse 7-ray emission, see [331 [3lj. 

3.2.1. Diffuse Galactic gamma rays. The diffuse emission (DE) from the Milky Way 
dominates the 7-ray sky in the LAT energy range. About 80% of the high-energy lumi- 
nosity of the Milky Way comes from processes in the interstellar medium (ISM). The 
Galactic diffuse 7-ray emission is the product of CR particle interactions with gas in the 
ISM and the low energy photons in the interstellar radiation field (ISRF). Therefore, its 
calculation requires first the estimation of CR spectra throughout the entire Galaxy |31j . 
Major components of the DE are vr^-decay, inverse Compton (IC), and bremsstrahlung. 

The first detailed analysis of the diffuse emission from the Galactic plane |6| < 10° was 
made by Hunter et al [35]. The spectrum of 7-rays calculated under the assumption that 
the proton and electron spectra in the Galaxy resemble those measured locally reveals an 
excess at > 1 GeV in the EGRET spectrum, the so-called " GeV excess" . An extensive 
study of the Galactic diffuse 7-ray emission in the context of cosmic ray propagation 
models has been carried out by Strong et al [311 [30]. This study confirmed that models 
based on locally measured electron and nucleon spectra and synchrotron constraints are 
consistent with 7-ray measurements in the 30 MeV - 500 MeV range, but outside this 
range excesses are apparent. In this paper the GALPROP "conventional" model will be 
used as a representative of this class of models. The resulting 7-ray spectrum from this 
model is shown in figure [31 

The "GeV excess" is seen in all directions, not only in the Galactic plane. A simple 
re-scaling of the components (vr^-decay, IC) does not improve the fit in any region, since 
the observed peak is at an energy higher than the vr^-peak. Although a possibility that 
the excess is an instrumental artifact due to the uncertainty in calibration can not be 
completely excluded [Ml 137], the detailed simulations by Baughman et al [38] indicate 
that the GeV excess is significantly larger when instrumental effects previously unac- 
counted for are considered. 

Alternatively, the assumption that locally measured electron and proton spectra 
represent the Galactic average can be dropped. In that case, as we will discuss in detail 
in the next paragraph, the "GeV excess" can be reproduced. The resulting model. 
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Figure 3. 7-ray spectrum of the inner Galaxy (300° < I < 30°,|6| < 5) derived 
from the "conventional" model (see text for more details). Dotted: contribution 
from TT*^ decay, dashed: contribution from inverse Compton scattering, dash-dotted: 
contribution from bremsstrahlung, solid: extragalactic background, bold solid: total 
flux. Also shown are the data points from EGRET (dark bars, red in colored versions) 
and COMPTEL (light bars, green in colored versions). Figure taken from [33] . 

dubbed "optimized GALPROP" model, is used as another benchmark model to represent 
possible diffuse 7-ray background. The resulting 7-ray spectrum from the "optimized" 
model is shown in figure HI 

The "optimized" model. There are a number of reasons why the CR intensity may 
fluctuate in space and time. First are the stochastic spatial and temporal distributions 
of supernova (SN) events. Dramatic increases in CR intensity, perhaps connected with 
nearby SN explosions, are recorded in terrestrial concentrations of cosmogenic isotopes. 
Concentrations of ^°Be in Antarctic and Greenland ice core samples [39] and ^°Fe in a 
deep-sea ferromanganese crust [ID] indicate highly significant increases of CR intensity 
~40 kyr and 2.8 Myr ago. The SN rate is larger in the spiral arms [H]. This may lead 
to lower CR intensity in the interarm region where the solar system is located. In the 
case of anisotropic diffusion or convection by the Galactic wind such fluctuations may 
be even stronger (see [31] for a discussion). The intensity variations and spectra of CR 
protons and heavier CR nuclei may be uncorrelated [12]. The total inelastic cross section 
for protons is ~30 mb vs. ~300 mb for Carbon, so that their Galactic "collecting areas" 
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differ by a factor of 10 or more for heavier nuclei. This means that carbon comes from 
Galactic regions that are about 10 times closer than the sources of protons, and thus 
implies that the directly measured CR protons and CR nuclei may come from different 
sources. 

Additionally, CR electrons and positrons suffer relatively large energy losses [43] as they 
propagate away from their sources and thus their spectral and intensity fluctuations can 
be considerably larger than those of protons. 

Antiprotons in CR are presumably secondary and produced mostly by CR protons in 
interactions with the interstellar gas. Because of their secondary origin in the ISM, 
which is far more uniform in space and time than the primary CR sources, their inten- 
sity fluctuates less than that of protons. The total cross-section of antiprotons is about 
the same as of protons, except at low energies due to annihilation, and thus they trace 
the CR proton spectrum averaged over a large region of the Galaxy. Therefore, if the 
directly measured local CR spectrum is not representative of the local Galactic average, 
then the antiproton measurements can still be used instead to derive this average inten- 
sity of CR protons. 

In addition to spatial diffusion in the interstellar space, the scattering of CR particles 
on randomly moving MHD waves leads to stochastic acceleration (second order Fermi 
acceleration or "reacceleration"), which is described in the transport equation as 
diffusion in momentum space. It has been shown [H] that the distributed acceleration 
may be strong enough to explain the peaks in the secondary-to-primary nuclei ratios at 
approximately 1 GeV/nucleon. When normalized to the local CR proton spectrum, the 
reacceleration model underproduces p's (BESS data [15]) at ~2 GeV and diffuse 7-rays 
above 1 GeV by the same factor of ~2 [30l US] while it works well for other CR nuclei. It 
is thus enough to renormaliz^ the CR proton spectrum up by a factor of 1.8 to remove 
the excesses. The model then predicts a factor of 2 too many photons at ~100 MeV. The 
100 MeV photons are produced mostly by ~1 GeV protons, where many uncertainties 
simultaneously come into play: poor knowledge of the 7r°-pro duct ion cross-section at 
low energies and/or low-energy interstellar proton spectrum and/or solar modulation 
(see jinillHllin] for further discussion). To get agreement with the EGRET photon data 
effectively requires a corresponding adjustment of the spectrum of CR protons at low 
energies. Since the IC and vr^-decay photons have different distributions, the electron 
spectrum also needs to be renormalized up by a factor of 4 in order to reproduce the 
EGRET diffuse 7-ray flux itself. These adjustments in the injection spectra of protons 
and electrons are enough to reproduce the spectrum of the EGRET diffuse 7-rays in all 
directions as well as the latitude and longitude profiles for the whole EGRET energy 

|j The BESS-Polar flight of 2004 [47] revealed that the CR antiproton flux is somewhat lower than 
previous measurements with lower statistics. This reduces allowable variations in the CR nuclcon 
spectrum due to the uncertainties in the antiproton flux measurements. 
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Figure 4. 7-ray spectrum of the inner Galaxy (300° < I < 30°,|5| < 5) derived from 
the "optimized" model (see text for more details). Dotted: contribution from 7r° decay, 
dashed: contribution from inverse Compton scattering, dash-dotted: contribution from 
bremsstrahlung, solid: extragalactic background, bold solid: total flux. Also shown 
are the data points from EGRET (dark bars, red in colored versions) and COMPTEL 
(light bars, green in colored versions). Figure taken from [33] . 



range 30 MeV - 50 GeV. 



3.2.2. Extragalactic diffuse •y-rays. The extragalactic 7-ray background (EGRB) is a 
weak component that is difficult to disentangle from the intense Galactic foreground. Its 
spectrum depends on the model adopted for the Galactic diffuse emission, which is itself 
not ffimly determined. The isotropic, presumably extragalactic, component of the dif- 
fuse 7-ray flux was first discovered by the SAS-2 satellite [50] and conffimed by EGRET 
[5T] . Even at the Galactic poles, the EGRB does not dominate, and its flux is com- 
parable to the Galactic contribution from IC scattering of the Galactic plane photons 
(interstellar radiation field) and cosmic microwave background photons by CR electrons 
[521 [53]. Additionally, Compton-scattered solar photons by CR electrons make broadly 
distributed emission with maximum brightness in the direction of the Sun [HU [551 156] . 
This all-sky average ~ 1 ■ 10~^ photons cm~^ s^^ sr^^ above 100 MeV is about 10% of 
the extragalactic emission (as inferred from EGRET data). Finally, the 7-ray albedo 
of small solar system bodies (asteroids) may contribute at high Galactic latitudes as 
an additional component of the 7-ray foreground distributed around the ecliptic [57] . 
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The determination of the EGRB is thus model-dependent an influenced by the assumed 
size of the Galactic halo, the electron spectrum there, and the spectrum of low-energy 
background photons, each of which must be derived independently. It can also be im- 
pacted by instrumental backgrounds, which are at a comparable level, in case they are 
not accurately determined and subtracted. 

More recently. Strong et al |52] reanalyzed the EGRET data using the "optimized" 
model discussed in section 13.2.11 Since this "optimized" model is not exact, the same 
method as in [51] was used. Reasonable agreement was obtained, considering the 
EGRET systematic uncertainty of ~15%. A larger deviation, was found only in one 
region, 2-4 GeV. This new estimate of the EGRB spectrum is lower and steeper than 
found by [51] in most energy ranges. It is not consistent with a power-law and shows 
some positive curvature, as expected in blazar population studies and various scenarios 
of cosmological neutralino annihilation [581 EH]. The integral flux above 100 MeV is 
(1.11 ± 0.01) ■ 10~^ photons cm-2 s"^ sr"^ [19]. 

The extragalactic background models, that are used as benchmarks in the present work 
are those presented by [51] and [52], which are consistent with EGRET data. The EBGR 
is most important in the analysis of cosmological WIMP annihilation (see section . 
For this analysis we also consider a model which takes into account GLAST's ability to 
resolve many of the sources that contributed to the EGRET EGRB measurement 

3.2.3. Particle backgrounds. Charged particles, mainly protons, electrons and 
positrons, as well as a smaller number of neutrons and Earth albedo photons, present 
a major instrumental background to potential DM signals. These background particles 
greatly dominate the flux of cosmic photons incident on the LAT, but the background 
rejection capability of the instrument is such that less than 1 in 10^ end up in the final 
sample of 7-ray candidates. The necessary cuts are included in the simulated instrument 
performance. 

We take into account the remaining charged particle background in the analyses where 
it can have a significant effect (i.e. for cosmological WIMP annihilation and the search 
for DM sateUites). 

4. Sensitivity to generic WIMP annihilation 
4.1. The Galactic Center 

The Galactic center (GC) is expected to be the strongest source of 7-rays from DM 
annihilation, due to its coincidence with the cusped part of the DM halo density profile, 
see Equation ([2]). As a result, the GC is frequently proposed for DM searches [16] and 
literature devoted to possible signatures of a DM 7-ray signal at the GC is extensive. 
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The "GeV excess" discussed in the previous section has been interpreted as the pos- 
sible signature of WIMP annihilation [231 [60], although such an interpretation relies 
heavily on the underlying model for the Galactic diffuse background (see discussion in 
section 13.21) . 

The neighborhood of the GC harbors numerous objects capable of accelerating CR to 
very high energies and thus producing 7-rays by inverse Compton scattering on elec- 
trons or pion decays following proton-proton or proton-7 interactions. Most notably, 
the bright, very high energy 7-ray point source observed by the H.E.S.S. [61], MAGIC 
[62], VERITAS [63] and CANGAROO ^3] collaborations represents a formidable back- 
ground for DM studies in this region of the sky. Initially proposed as the DM signature 
of a very heavy (> 6 TeV) WIMP ([6^ and references therein), such a signal would rule 
out a Moore profile for neutralino masses above ~400 GeV, for < cxw > ~ 10~^^ cm'^ s~^ 
or the expected cross-section for Kaluza-Klein particles [66]. However, this source is 
now widely considered to be a standard astrophysical source, either associated with the 
~ 2.6 ■ IO^Mq black hole at the kinematic centre of our Galaxy, commonly identified 
with the bright compact radio source Sgr A*, or with the candidate pulsar wind nebula 
G359.95-0.04 [67] recently discovered in a deep Chandra survey and only at 8.7" from 
Sgr A*. One should also take into account the presence of a possible EGRET source at 
approximately 0.2° from the GC [68] . 

At this point, we do not know the behavior of these high energy sources in the GLAST 
energy range, which makes extrapolation for the sake of including the source flux as an 
additional background doubtful. The following analysis is based on the hypothesis that 
the main task will be to distinguish the DM signal from a Galactic diffuse background 
after the astrophysical sources are disentangled and subtracted using the information 
provided by spectral and angular analysis and multiwavelength observations. This sub- 
traction will introduce a systematic uncertainty in the assessment of the diffuse flux 
which we neglect here. For a recent attempt to include HESS and EGRET sources into 
a sensitivity calculation, see [69] . 

As described in section [3l the expression of the 7-ray continuum flux for a generic WIMP 
at a given energy E is given by Equation ([1]). A truncated NEW profile as defined in 
[23] is assumed for the WIMP distribution, and only one dominant annihilation channel 
(ly+iy-^ 66, tt, r+r^) is considered at a time. Care has also been taken in order not to 
violate the EGRET flux constraint around the GC [70]. WIMP annihilation differential 
fluxes above 1 GeV and in a region of 0.5 degs radius around the GC (corresponding 
to the angle for 68% containment at this energy threshold) have been generated using 
DarkSusy v. 4.15 [71]. The expected DM from this region, incident on the LAT for 5 
years of all sky scanning operation is simulated using gtobssim (see section [2^ . This 
flux was simulated for 10000 grid points (mwiMP, < o'v >) on a 100 x 100 logarithmic 
grid. The range of chosen is < av >e [10^^^, 10"^^] and mwiMP e [10, 3000] GeV for the 
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bb and r+r~ channels, and mwiMP ^ [200, 3000] for the W^W~ and ti channels. The 
EGRET flux modeled in [70] is also simulated at the GC and a standard statistical 
analysis is performed to check if a given WIMP model conflicts with EGRET data at 
the 5(7 level. For models compatible with EGRET data, a second test is performed 
to check if GLAST is able to disentangle the WIMP contribution from either the "con- 
ventional" or the "optimized" Galactic diffuse background. For the calculation, 20 
logarithmic bins in the energy range [1 GeV, 100 GeV] are used for hb and r+r^, [1 GeV, 
300 GeV] for W^W~ and ti. Figures [5] to [7] show the results of the scans for different 
annihilation channels. These results are given at a 3(j confidence level. We also show 
the 5cr significance in figure [5] for the bb final state. As expected, sensitivity regions for 
the 66, ti and W^W~ annihilation channels are quite similar, whereas differences are 
apparent for the r+r^ annihilation channel^. Furthermore, it can be seen that regions 
of O(10'26)pj^3g-i g^^g within the reach of GLAST. 

In the GC study presented here, a NEW profile is assumed. Recently, Hooper et al 
[72] have studied the observational consequences for GLAST of an excess in the WMAP 
foreground coined the "WMAP haze" [73] , when interpreted as the synchrotron emission 
from highly relativistic electron-positron pairs produced by WIMP annihilation. The 
author's conclusions are encouraging. If a slightly steeper inner slope of the NEW profile 
is considered, p(r) oc r~^'^ rather than r~^, as suggested in [68], the LAT sensitivity 
should be improved by a factor of ~ 10. Assuming a Moore profile, the line of sight 
integral increases by a factor of 200 (with corresponding improvement in sensitivity), 
under the assumption of cored profile the sensitivity decreases by a factor of ~100. 
Compared to this, a simplified introduction of a 20 % uncertainty in the Galactic diffuse 
background (by adding a term to the x^), only leads to a decrease in sensitivity by about 
5 %. It is noteworthy, that observations of the GC at x-ray and radio wavelengths might 
lead to strong constraints on the DM signal [71] . 

4-2. The Galactic Halo 

The detection of the continuum WIMP signal from the Galactic halo is complicated by 
difficult systematic effects due to the presence of comparably uncertain astrophysical 
backgrounds (see section [372]) . However, analyzing the halo is complementary to analy- 
ses of the immediate region of the GC, where the DM rates are significantly greater in 
the standard halo models, but where there are also issues of source crowding and large 
backgrounds. A potential continuum WIMP signal has the advantage that its spectral 
shape would be the same no matter where observed in the Galaxy assuming effects of the 
astrophysical environment can be neglected. In addition, in most scenarios one expects 
a sharp cutoff in this signal at the mass of the WIMP that is very difficult to achieve 
with expected astrophysical backgrounds. 

^ For the other analyses presented in this paper, we usuaUy assume hh as the dominant annihilation 
channel. Strictly, those regions might not be useful for models with dominant annihilation to tt^ . 
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Figure 5. Cross-sections < av > {v ^ 0) versus the WIMP mass mwiMP for the 66 
annihilation channel. Left panel shows the result for 3a significance, right panel shows 
the result for 5a significance for 5 years of GLAST operation, The upper part 
of the plots corresponds to regions which are already excluded by the EGRET data 
around the GC and the lower part corresponds to regions not detectable by GLAST. 
The "detectable by GLAST region" corresponds to models detectable by GLAST for 
both "conventional" and "optimized" astrophysical background. The shaded region 
represents models which can be detected only under the assumption of "conventional" 
Galactic diffuse background. See text for more details. 



This section presents an analysis method for the WIMP continuum contribution to the 
Galactic diffuse emission in the halo. For this analysis we include the full Galactic Halo 
except for a region excluded because of large diffuse backgrounds and systematic errors 
associated with this background. We consider two possible exclusions: the first is the 
region within 10° of the GC and the second is the region within 10° of the Galactic 
plane. Both of these choices will avoid the GC where the proposed dark matter density 
profiles diverge from one another. 

The statistical sensitivity to a WIMP signal is quantified as a function of the WIMP 
mass and other parameters such as the annihilation cross-section and DM density profile. 

Since the diffuse 7-ray emission is the largest background for a halo DM search, two 
models ("optimized" and "conventional") of the Galactic diffuse background intensity 
maps (see section 3.2) are used in the analysis. We generate DM models using a NFW 
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Figure 6. 3 cr sensitivity regions for tlie tt (left panel) and the Vt^^ annihilation 
channel (right panel) . Definitions of regions are the same as in figure [51 Note the 
difference in x-axis scale as compared to figure [5l 



density profile. With the intensity maps calculated from the diffuse model plus the DM 
model computer simulations of the detector are performed to obtain GLAST events with 
the expected distributions in space, time and energy. Figure |8] shows the shapes of the 
energy distributions observed in the LAT for various DM masses versus the two models 
of the diffuse background, "conventional" and "optimized" . 

The Galactic diffuse background peaks at about 150 MeV, while the DM distribution 
peaks at about 1 GeV for a WIMP mass of 50 GeV/c^, and increases with increasing 
WIMP mass. Background and signal spectra are thus distinguishable. In addition, the 
DM has a different spatial distribution. We can take advantage of these differences by 
performing a simultaneous fit over angles (i.e. Galactic latitude and longitude) and in 
energy. We perform a log-likelihood fit of the dark matter and diffuse model to nu- 
merous pseudoexperiments generated from the same models binned 3 x 3 x 30 in /, 6, 
and logE). In the fit there are no constraints to the normalization of either signal or 
background. 

We estimate the LAT sensitivity for an observation of DM from WIMP annihilation 
by running sample experiments (or pseudo-experiments), assuming data sample sizes 
consistent with one year of running. Each pseudo-experiment includes photons from 
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Figure 7.3 a sensitivity regions for the t+ t annihilation channel. Definitions of 
regions are the same as in figure [S] 
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Figure 8. The energy distributions for the diffuse "conventional" and "optimized" 
diffuse-emission models and two DM models for WIMP masses at 50 and 250 GeV/c^. 
Each distribution is individually normalized to unit area to highlight the shape 
differences. 
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WIMP annihilation with fixed mass and an assumed < av >, along with photons ex- 
pected from diffuse gamma ray background. Other backgrounds are smaller and are 
presently ignored. We then fit the resulting pseudo-experiment to the background and 
signal shape in latitude, longitude, and energy. We construct ensembles of ~ 1000 such 
pseudo-experiments, with each ensemble corresponding to a fixed WIMP mass from 50 
to 250 GeV/c^. For each mass, we vary < av > until 50% of the experiments have 
observations at either 3a or 5a significance. Figure [9] shows the results for both the 
"conventional" and the "optimized" diffuse backgrounds. 

In our fiducial region, a Moore profile predicts a factor 1.3 times the fiux of NFW, which 
would enter linearly into our sensitivity. While the uncertainty in the signal fiux (due to 
lack of knowledge on DM density distribution) is the dominant systematic uncertainty 
entering the sensitivity calculation, the question if a halo DM signal is distinguishable 
from a the cosmic-ray induced 7-ray background will crucially depend on the background 
uncertainty, for which careful studies with data will have to be performed and which 
requires sophisticated statistical tools due to non-linear dependences in the parameters 
entering the background model. 

Nevertheless, in order to estimate the impact of a systematic uncertainty on the shape 
of the diffuse background in a rough fashion, we introduce a nuisance parameter into the 
likelihood function. The parameter allows for a continuous change in the background 
distribution from the "conventional" model to the "optimized" model. The likelihood 
is minimized for the nuisance parameter and the signal and background components. 
A fiat prior is assumed for the nuisance parameter. When using the "conventional" 
background model and allowing for a systematic variation in shape, primarily in the 
energy spectra, the required < av > for a 3a {5a) observation is increased by ~70% 
(~45%) at low mass and ~40% (~25%) at high mass. When using the "optimized" 
background model the change due to the systematic variation of shape is less, about a 
~30% change at low mass and a ~20% at high mass for both the 3a and 5a observations. 

It can be concluded, that (with the caveat of only preliminary treatment of systematic 
uncertainties) the LAT would probe a large region of the MSSM and mSUGRA param- 
eter space, as shown in figure [25] in section [51 

Using a similar technique we obtain the expected LAT mass resolution provided that 
< av > would be at a high enough level. Again, we construct ensembles of 1000 pseudo- 
experiments, this time with fixed WIMP masses and < av >= 2.3 ■ 10~^^cm^s~^. For 
an observation time of 1 year, this results in approximately 3.3 -10^ photons from a 250 
GeV WIMP signal, and 10^ photons from diffuse background ("conventional" model). A 
likelihood fit of each pseudo-experiment is made to a combination of diffuse background 
and signal, varying the assumed WIMP mass. An example of such a fit for a single 
pseudo-experiment is shown in the left side of figure [TOl The resulting fitted mass for 
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Figure 9. The < av > required to obtain an observation of WIMP annihilation at 
either Sa (square) and 5cr (circle) significance for one year of GLAST data, as a 
function of WIMP mass. Left panel: considering the "optimized" diffuse model as 
background, right panel: considering the "conventional" diffuse model. The dashed 
line corresponds to the 10° cut above and below the Galactic plane; the solid line 
corresponds to a 10° radial cut around the Galactic center. 



this pseudo-experiment is 155±23 GeV/c^. We then repeat this procedure for additional 
WIMP masses of 100, 200 and 250 GeV/c^. Assuming the same av = 2.3 ■ 10"2^cm^s"^ 
and 1 year of observation resuhs in signal counts which vary from 1.5 ■ 10^ (at 100 
GeV/c^) down to 3 ■ 10^ (at 250 GeV/c^). The right side of figure [TU] shows the average 
returned uncertainty on the fitted mass versus WIMP mass. From this figure, we see 
that our expected mass resolution varies from about 11 GeV/c^ (at 100 GeV/c^) to 63 
GeV/c^ (at 250 GeV/c^). The bars in figure [TOl show the 68% containment interval 
for possible outcomes of a single experiment. The mass resolution shows that from the 
consideration of the statistical errors it would be possible to resolve the mass of the 
WIMP, which could be used to guide further 7-line and accelerator searches. Note that 
the mass resolution obtained in this analysis is not dominated by the energy resolution 
of the detector (in contrast to the case of detection of a line signal), but rather intrinsic 
to the signal. 

4.3. Galactic Satellites 

4.3.1. Galactic sub-halos. In the CDM paradigm [75] [76], structure forms hierarchi- 
cally, so the dark halos of galaxies such as the Milky Way should contain large numbers 
of subhalos. For WIMPs the subhalo mass spectrum is expected to extend as low as 
10~^ to 10~^^ solar masses depending on the free-streaming scale of the considered DM 
particle, see [77]. It is improbable that GLAST will detect halos of this size as they 
would need to be very close to the Earth. 

The substructure is expected to be nearly isotropic. Thus, annihilation in the subhalos 
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Figure 10. The left plot shows the fitted —A log-likelihood versus mass for a 
single pseudo-experiment assuming the true mass of the WIMP to be 150 GeV/c^. 
The resulting fitted mass for this pseudo-experiment is 155 ± 23 GeV/c^. The right 
plot shows the mean fitted error returned by the fit, versus the true WIMP mass. 
The error bars show the 68% containment interval for the pseudoexperiments, which 
demonstrates the range of possible outcomes from a single experiment. The parameters 
describing the pseudoexperiments are given in the text. 



can be significant away from tlie Galactic plane, where astrophysical sources are con- 
centrated. These substructures are commonly called DM satellites. 

To take into account the tidal stripping of DM satellites, we use a truncated NFW 
profile [78] , which approximates a structure of a stripped halo by a simple modification 
of the NFW profile: 



p(r) 



where r is the distance from the center of the structure, ft is a dimensionless measure of 
the reduction in central density, and is an "effective" tidal radius. A Milky Way-sized 
dark halo is simulated based on the above satellite profile and a satellite distribution 
generated by Taylor & Babul [79], which was chosen as representative among many 
possible theoretical predictions (see for example [80], [81]). In their simulations, about 
30 % of the mass is in satellites. Then, a generic WIMP model, in which the WIMP 
mass is 100 GeV, the annihilation cross-section is 2.3 ■ 10~^^cm^s~^ (consistent with 
the WMAP DM density estimate) [82], and the annihilation channel is 66, is used to 
estimate the number of Milky Way DM satellites with > 10^ solar masses observable by 
the LAT after 1 and 5 year exposures in scanning mode, described in section 2.1. The 
DM satellite distribution is roughly spherically symmetric about the GC and extends 
well beyond the solar orbit. Thus, the observable satellites are located mostly at high 
Galactic latitudes. The background is estimated using the isotropic extragalactic dif- 
fuse background described in Sreekumar et al [5lj plus the GALPROP "conventional" / 
"optimized" model as discussed in section [3^ The significance of the DM signal is then 
estimated to be the number of signal events within the satellite tidal radius (or the PSF 
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68% containment radius, whichever is larger) divided by the square root of the number 
of total events, including both background and signal, within the same radius, at energy 
larger than 1 GeV. Figure [TT] shows the cumulative number of DM satellites with sig- 
nificance of at least 5 a. For this generic WIMP model, the LAT is expected to observe 
~12 such DM satellites within 5 years, (assuming uniform exposure). For discussion 
on observability for other models see for example [83], [H]. Adding the charged parti- 
cle background to this analysis, reduces the number of detectable satellites by about 
10 %. Uncertainties in the charged particle background prediction have negligible effect. 

We also calculate the LAT error ellipses for WIMP mass vs. the quantity given in 
Equation (jlj) for this generic WIMP model. We choose one "5 a" satellite for 1 year 
of GLAST data (this significance is estimated using the back of the envelope method 
mentioned above), which is at high latitude (b = -39 deg), 8.9 kpc away from the Earth, 
with the mass 2.7- 10'' M©, and with the tidal radius 0.2 kpc. We use the likelihood ratio 
test statistic and profile likelihood to extract the 99%, 90% and 68 % error contours on 
WIMP mass and annihilation cross-section jointly, as shown in figure [121 fo^^ this generic 
WIMP DM satellite, where the unit of the vertical axis is given by 

<av> /8.9kpcV / MsatcUitc \^ /0.2kpc^^ 



2.3 X 10-26 cm3s-i ^ V d ) ^ \2.7 xWMq) ^ V / ^^"^ 

which can be used to rescale the error ellipse for different satellite masses, distances 
and tidal radii. The above scaling factors are calculated from the gamma ray flux from 
WIMP annihilation, for point-like sources, for which 

<PwiMP oc Vf{V) (5) 

X 

where 



is a dimensionless flux multiplier [85] and p is the average density within volume V. 
In Taylor and Babul's simulations, the observable satellites for the LAT have ^ ~ 0.1, 
therefore it is accurate enough to consider f(V) over the volume within rt. For the NEW 
profile, for which p oc in that region, f(V) is a constant. Consequently: 

< av > 1 <aV > 1 Msatellite /-x 

'^^^ ^ -^d-^p ^ = -^d^^^ 

However, the Moore profile p7J with the same scale radius and scale density, exhibits 
p oc r~^-^ within r^. In this case, fiV) oc In rj/rmm- Integrating Equation ([5]) from r^j^ 
to Tt, taking into account the change in p, yields that the example satellite will emit ~50 
times more photons under the assumption of a Moore profile than under the assumption 
of a NEW profile (for a cutoff radius of r^m = 10"^ kpc). It should be noted that the 
Moore profile is rather optimistic. 
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In order to identify DM satellites, we need firstly to find possible candidates for DM 
satellites; secondly to make sure these candidates are not fluctuations of the background 
emission; and finally to distinguish these real candidates from other typical astrophysical 
sources, such as pulsars or molecular clouds. The typical extent of the observable DM 
satellites is on the order of 1 deg. For bright enough sources, the LAT should be 
able to resolve them using WIMP annihilation photons above 1 GeV, for which the 
PSF is < 0.5 deg. SExtractor [86], which can build a catalog of objects from an 
astronomical image, is used to search for possible candidates for DM satellites. We 
then use the binned likelihood analysis (see section 12. 2p to reject false sources and other 
typical classes of sources. Figure [13] shows the counts spectra for the same generic 
WIMP satellite, plus diffuse background. The diffuse background model consists of the 
GALPROP "conventional" / "optimized" model and the isotropic extragalactic diffuse 
background [51j • In future data analysis we will use a hypothesis test using a likelihood 
ratio test statistic to distinguish the satellite from background. As a sanity check of 
the methodology developed by the GLAST collaboration, we generate 120 simulated 
GLAST experiments for the background only and also for the example satellite. Then, 
we compute our test statistic (TS) as the maximum log-likelihood difference between 
background only and background plus signal, where the free fitting parameters are 
the signal and background normalizations. Figure [H] shows a two dimensional plot 
of the TS for background only (filled triangles) and for background plus signal (filled 
squares) versus best-fit signal strength, which is the ratio of the fitted flux to the actual 
simulated flux of the satellite. Indeed the TS seems to behave as expected. The signal 
satellite can be detected with an efficiency close to 100 %. Only one (two) out of the 
120 realizations gives a TS value which falls below the threshold of TS = 25 for the 
"conventional" ("optimized") GALPROP background. Using the unique WIMP pair 
annihilation spectrum, which is an extremely hard non-power-law and has an end point 
at the mass of the WIMP, we will use the same TS to distinguish DM satellites from 
other astrophysical sources. Also, note that real DM satellite candidates are expected 
to have no counterparts in radio, optical. X-rays and TeV energies. For more details on 
separating DM satellite signals from possible backgrounds, see [87] . 

4-3.2. Dwarf Galaxies. As described in the previous section, N-body simulations pre- 
dict a large amount of substructure which represent good targets for DM annihilation 
7-ray studies with the LAT. Dwarf spheroidal galaxies (dSph) may be manifestations 
of the largest clumps predicted by the CDM scenario. They are potential candidates 
for indirect searches of DM because they are apparently among the most extremely DM 
dominated environments. For example, the mass-to-light ratio of Draco is ~ 250 in solar 
units [90], while it is ~ 100 [91] for the Sagittarius Dwarf. 

As an illustration of the LAT capability to detect DM 7-ray signals from dSph galaxies, 
we compute the 3a sensitivity for 5 years operation for Sagittarius, which is the closest 
dSph to the Sun. Discovered in 1994, Sagittarius is one of the nine dwarf spheroidal 
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Figure 11. Estimated number of observable DM satellites for the LAT in the 
Milky Way for 1 and 5 year of GLAST operation. The background consists 
of the isotropic extragalactic diffuse [5T] and GALPROP "conventional" / "optimized" 
Galactic diffuse model. The significance is estimated as Ns / V Ns + Nb within the 
satellite tidal radius (or the PSF 68% containment radius if larger ) at E'^.y > 1 GeV. 
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Figure 12. LAT error ellipses for a simulated "5 sigma" DM satellite with the generic 
WIMP model. The vertical axis is given by Equation ([4]). 



galaxies orbiting our own Galaxy. The core of the Sagittarius dwarf is located at 1=5.6" 
and b=— 14'' in Galactic coordinates at a distance of about 24 kpc from the Sun. As 
most dSphs do not contain significant amounts of gas, the structure of the dark halos 
must be inferred from stellar motion. For this study we use the Moore cusped halo 
profile derived for Sagittarius in [SB]. For compatibility with EGRET observations, we 
require here that the total integrated flux given in Equation ([1]) does not exceed the 
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Figure 13. Counts spectra for the same generic WIMP satellite within its tidal radius. 
The squares with 1 sigma error bars show the spectrum of simulated 1 year GLAST 
data. The thick lines show the background plus the satellite signal. The thin lines 
show the satellite predictions only. The dash-dot lines show the background predictions 
only. 
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Figure 14. Null hypothesis test to distinguish the same generic WIMP satellite from 
a background fluctuation. The test statistics for background only and for background 
plus satellite signal could be separable. 



value(0.4 ■ 10 cm ^ s ^) extracted from the upper limit maps presented in \89\ . 

The resuh reported in figure [I5] shows the LAT capabihty to detect a DM signal from 
the Sagittarius dSph in the case of the favorable cusped Moore profile proposed in [88] . 
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Figure 15. Sensitivity to a Sagittarius Dwarf DM signal for 5 years of GLAST 
operation assuming Moore profile as described in [88]. The region labeled "above 
EGRET observation" is calculated with respect to the upper limit map provided in 



Using a Moore profile for Sagittarius is a very optimistic assumption, though it has 
been argued that the number of passes through the Galactic disk is low enough, so that 
it is not completely disrupted [88]. More conservatively assuming a NFW profile will 
decrease the sensitivity by about a factor 10. 

For example, the assumption of a Moore profile is less optimistic for the Draco dwarf. 
However, Draco is about a factor 3-4 more distant. Correctly taking into account the 
difference in Galactic diffuse background, leads to an estimated decrease in sensitivity 
by about a factor 10, when comparing Draco with Sagittarius assuming a Moore profile. 
Since Draco on the other hand might have a more favourable densitiy profile, it is one 
of the prime candidates for DM searches. 

In conclusion, nearby dwarf spheroidals are attractive targets for DM searches with the 
LAT observatory, if cuspy profiles can be considered. 

4.4- Point Sources 

In this section we study the possibility for GLAST to detect point sources of Dark 
Matter. These sources can be motivated by scenarios invoking intermediate size black 
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holes [22], i-e. wandering black holes with masses 10^ < M/Mq < 10^, which would 
adiabatically grow "mini-spikes" of DM. 

In order to study the LAT sensitivity to such objects, we perform a two- month scanning 
mode simulation of the 7-ray sky. For the backgrounds we employed the "optimized" 
GALPROP model for Galactic diffuse emission and the EGRET measurement [51] for 
the extragalactic contribution. We divide the sky into 23 regions of about 10 degs in 
radius, and in each region we placed one DM point source. Then, we consider each 
source separately and let the flux intensity above 100 MeV vary from 10^*^ to 10"'' ph 
cm~^s"^. These flux values are typical for scenarios with mini-spikes as for example 
discussed in [92]. For each intensity, we calculate the significance of the observed signal, 
given the local background counts, with a maximum likelihood analysis assuming Pois- 
son statistics as described in section 12. 2[ By estimating the minimum flux required to 
discriminate the DM source from the background at a 5(j level on a grid of points uni- 
formly distributed over the sky, we obtain the sensitivity map shown in figure [16] (where 
we adopt a DM particle mass my/iMP = 150 GeV and assume annihilation mainly into 
hi). 

The sensitivity appears to depend significantly on the Galactic longitude only along the 
Galactic disk, as expected. At high Galactic latitudes a source as faint as 1 ■ 10~^ ph 
cm~^s~^ above 100 MeV is resolved, while close to the GC a minimum flux of 8 ■ 10~® ph 
cm~^s~^ is required. Figure [H] shows some illustrative examples of simulated sources. 
All sources detected with a significance above 5 a show a clear indication of a spectral 
cutoff at high energy, i.e. the fit with a simple power line spectrum is disfavored. Mini- 
spikes scenarios such as the one discussed in [92] predict a population of ~ 100 DM 
mini-spikes. Using the likelihood approach to discriminate against usual astrophysical 
background we expect to be able to unambiguously detect about 25 % of the sources 
(depending on mini-spike distribution and assumptions on background). 

4-5. GLAST sensitivity to a line signal from DM annihilation 

In this section an approach for the indirect astrophysical detection of DM annihilation 
lines is presented. A line at the WIMP mass, due to the 27 production channel, could 
be observed as a feature in the astrophysical source spectrum [16] . Such an observation 
is a "smoking gun" for WIMP DM as it is difficult to explain by a process other than 
WIMP annihilation or decajll]. In addition, it should be free of astrophysical uncertain- 
ties, since the background can be determined from the data itself. Remember however, 
that in the most popular models branching ratio for the annihilation into lines is typi- 
cally about 10"^ or less (see [17], however). 

We calculate sensitivities here in units of photons/cm^ /s/sr which thus are independent 

+ The presence of a feature due to annihilation into 7Z in addition would be even more convincing. 
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Figure 16. Flux needed for separating a point sources of DM anniliilation from the 



background, i.e. full-sky map in Galactic coordinates of the minimum flux above 100 
MeV, in units of [ph cm^^s^^], that is required for a 5tT detection of an annihilation 
spectrum, assuming a DM particle with mass tow/a/p — 150 GeV annihilating into 
bb (note, however, that the map does not depend very sensitively on DM properties). 
The map is relative to a 2 month operation period; for longer operation times, 
fluxes scale approximately as t^^g"^ ■ 



of the signal model. However, they do depend on the diffuse background model. The 
Galactic diffuse models used here are the GALPROP "optimized" and "conventional" 
models, discussed in section 13. 2[ An isotropic power law is used for the extragalactic 
diffuse component [51]. Theoretical predictions can be compared with these sensitivities 
assuming a WIMP model (e.g., WIMP mass and many other theoretical parameters), 
a DM halo model (e.g. the NFW profile) and other astrophysical backgrounds. To 
estimate the photon background, this analysis uses diffuse-model simulated data from 
the Galactic centered annulus(r G [20°, 35°]), excluding the region within 15° of the 
Galactic plane. In the context of an NFW DM profile [26], including data in such 
an annulus minimizes the contribution of the Galactic diffuse emission and could give 
a signal-to-noise ratio as much as 12 times greater than at the GC [93]. The region 
we use in this analysis is not necessarily optimal for real data, but is shown as an 
example of how one might optimize signal to background. LAT line energy sensitivities 
are calculated for 5cr detection and for the 95% confidence level upper limit (CLUL). 
We address cases for which the line energy is known (e.g. supplied by a discovery at 
the Large Hadron Collider [LHC]), as well as for which the line energy is unknown, 
requiring a search over the energy range of interest. The latter gives somewhat poorer 
limits as the number of energy bins introduces a trial factor. The LAT resolved signal is 
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Figure 17. Examples of spectral fits of simulated DM point sources of intensity 
(f> = 8 ■ 10"^ photons cm^^s'^ above 100 MeV, miy/MP = 150 GeV 66 annihilation 
channel. Upper left for (l,b) = (0,20); upper right for (l,b)= (0,50) and on the bottom 
(l,b)=(60,0). Thin solid lines: DM signal, dashed: Galactic diffuse contribution, 
dotted: extragalactic contribution (from |51j). points with error bars are photon counts 
from the simulated observation. 



simulated using approximate delta functions as input for the full detector Monte Carlo 
(MC) simulation and reconstruction framework, GLEAM and applying event selections 
corresponding to the ones developed by the LAT collaboration. To calculate known line 
energy sensitivities, LAT resolved lines are generated from 50 GeV to 300 GeV at 25 
GeV intervals. The results for uniform exposure over the LAT are reasonably well fit 
by a double Gaussian distribution 0i: 



i{E;Ei^2,cri,2,NT,r) 



cr2 



{E-E2f/2al 



(8) 



where Nt = Ni + N2, and r = N2/NT. We use a double Gaussian instead of the of- 
ficial parameterization of the energy dispersion, since it is simpler to handle and since 
it gives a good fit over the entire energy range of interest. The center of the signal 
Eo is defined by the peak position. Figure [18] shows some examples of fits. The full 
width at half maximum (FWHM) is used to define an equivalent single Gaussian aE for 
each line. We use gtobssim to generate LAT resolved 5 year all-sky diffuse background 
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Figure 18. Double Gaussian fits to tlie simulated LAT energy dispersion for 7-ray 
energies of a) 50 and b) 150 GeV. 
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Figure 19. "Conventional" Galactic plus extragalactic background fit to a powerlaw 
over the range [Eq ~ Gcte, Eq + 6<7e], for a) Eo = 75 GeV and b) Eq = 225 GeV. The 
fit range is shown as a shaded area. 



spectra in the region of interest, the Galactic annulus excluding the Galactic plane as 
defined above. Figure [12] shows, for Eq = 75 and 225 GeV, photon count spectra for 
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Figure 20. "Optimized" diffuse background and a 5a signal at 200 GeV. The black 
dots and open squares correspond to the diffuse background and the diffuse background 
plus MC signal, respectively. Full and dotted lines correspond to the signal plus 
background fit to 02 and 01+02, respectively. < Ax^ >=25.0 for this run. 

this background, together with an exponential fit {(j)2{E; a,b) = a ■ e^^^) over the range 
[Eq — Q(Te, Eq + 6(7e]. The spectra are well fit over this energy range. 

Next, the 5 year 5cr signal sensitivity is estimated. For each line energy the in- 
put background is resampled (bootstrapped) 1000 times with a 0i MC signal and 
fit successively to (0i + 02) and (/)2 with free parameters a, b and Nt, for the range 
[Eq — 6aE,EQ + 6aE (see figure [19]). Each bootstrap is randomly sampled from the 
original background resulting in mildly correlated background realizations. This series 
of 1000 bootstraps is rerun varying the number of events thrown into the 0i signal MC 
until < Ax^ >=< X^2+0i ~ '^)' example of a line signal just fulfilling 

this condition is shown in figure [201 

The average number of signal counts needed at each energy is then converted to the 
known line energy sensitivity using average exposures over the annulus (see figure [2T]) . 
The LAT detection sensitivity for a line of unknown energy in the range [40,350] GeV 
for a 5a above background signal corresponds to a confidence level of 99.99997%. To 
calculate the number of counts needed when the line energy is unknown, the probability 
of no signal detection in a single bin is used: 

1 1 

(1-P)"M- = ^= / e~ dx , (9) 
v27r J-oQ 



34 




Figure 21. 5a sensitivity contours (5 years of GLAST operation) in a) flux 
and b) velocity-averaged effective cross-section. Triangles and squares correspond to 
the "conventional" and "optimized" Galactic background model, respectively. Full and 
dotted lines correspond to the case of known and unknown WIMP energy, respectively. 



where P = 3 ■ 10"^ is the probabihty of detecting a signal greater than the "number 
of a", Her, in one or more of ribins energy bins. A bin width of the LAT FWHM/E = 
0.19, (cr/E = 0.08), was used, giving 15 bins over the energy range [40, 350] GeV. The 
value of the number of as needed for a 5a detection over 15 bins is found to be 5.5. The 
unknown line energy sensitivity is estimated by scaling the known line energy sensitivity 
by (5.5/5)2 (^ggg figure EH). 

To obtain the effective < av > in the WIMP annihilation, we assume an NEW profile 
and calculate the line of sight integral (see Equation ([T])) within the annulus region with 
the mass of the WIMP equal to the line energy in each bin that we considered. The 
right panel in figure EH shows the 5a cross-sections. 

The 5 year 95% confidence level upper limit sensitivity for known line energies is ob- 
tained similarly to the 5a case by bootstrapping the diffuse background (with no MC 
generated signal), and by fitting each bootstrap sample to 0i -|- 02 with free parameters 
a, b and Nt, for E G [Eq — 6aE,Eo + Gas]- Nt ± aNr is retrieved from the fit, and 
the 95% CLUL is obtained as l.GAajyT- The case of unknown hne energy {E G [40, 350] 
GeV) 95% CLUL is calculated as in the 5a case with scale factor 2.71/1.64. Over one 
thousand bootstraps were required to produce 1000 convergent fits for Eq = 50 and 75 
GeV. Results are shown in figure E21 
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Figure 22. 95% confidence level upper limit contours (5 years of GLAST 

operation) in a) flux and b) velocity-averaged efTective cross-section, as a function of 
the WIMP energy. The legend is identical to figure [211 



These results are preliminary as they lack the WIMP continuum contribution to 02, the 
event selections and corresponding parameterizations of the response are being updated, 
the optimum (/, b) range of integration is still being considered including the GC, and 
refined statistical methods utilizing Poisson statistics are being explored. Note that 
with the annulus and the NFW profile used in this paper, boosts of 500 or more are 
needed to set limits or see signals in currently interesting MSSM cross-section ranges. 
With the (optimistic) Moore profile there is not much change in the result for the region 
of the sky integrated over. This is because the lower boundary of this region is 15 ° 
from the GC, and the line of sight integral does not probe the cuspy region of Moore vs. 
NFW. Indeed, an isothermal-cored profile also give about the same result in this case. 
However, as already mentioned one should note that models with enhanced branching 
fraction into lines have recently been suggested jT7] . 



4-6. Cosmological WIMP annihilation 

WIMP pair annihilation, taking place in DM Halos at all redshifts might contribute 
to the extragalactic 7-ray background (EGRB, see section [3.2.2p . The 27-ray channel 
could then result in a distinctive feature, a line which is distorted by the integration 
over all redshifts. The signal of cosmological WIMP annihilation is subject to significant 
astrophysical uncertainties but less sensitive to the exact choice of halo profile due to 
the integration over large volumes. 
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There are several ingredients necessary to calculate the flux of 7-rays from cosmological 
WIMP annihilation. In addition to the 7-ray yield per annihilation (see Equation 
([1])), assumptions need to be made on the distribution and structure of DM halos on 
cosmological scales as well as on the model for the evolution of the Universe. One also 
needs to take into account absorption of the high-energy 7-rays on the extragalactic 
background light. Following [91], the flux of photons can be calculated as 

c?0^ < av > c pI fin ^s^'^i^) d.N^iEo{l + z)) 



J h[z) dE 



dEo 8n Hq rri^ 

where c is the speed of light, Hq the Hubble constant, < av > the annihila- 
tion cross- sect ion, mwiMP is the WIMP mass, po the average DM density, h{z) = 
A/nM^rH^^y3~+T^^^TH^^)M-TyA parameterizes the energy content of the Universe, 
A^(2;) describes the averaged squared over-density in halos as a function of redshift 
and t{z, Eq) is the optical depth. The 7-ray yield is the same as discussed in section 
I3.lt where we assume dominant annihilation into bb and a branching ratio of 5 ■ 10~^ 
into lines. 



The extragalactic 7-ray signal is strongly affected by absorption in the inter-Galactic 
medium, especially at high energies. The dominant contribution to the absorption in 
the tens of GeV-TeV range is pair production on the extragalactic background light 
emitted in the optical and infrared range. For the optical depth as a function of both 
redshift and observed energy we use the results of [95]. Newer results imply a slightly 
lower optical depth at low redshifts and a slightly higher depth at high redshifts, which 
in turn slightly enhance and weaken our signal at high and low energies respectively 
[96]. The cosmic evolution of the energy content of the Universe, h{z), is estimated 
from the WMAP three-year data Example spectra of a cosmological WIMP sig- 
nal together with the backgrounds considered in the present note are shown in figure [23l 

The structure of the DM density is encoded in the quantity A'^{z) in Equation ffTOj) . 
To derive this quantity one has to integrate the contribution from halos of all masses, 
weighted by the halo mass function, and for each mass take into account the spread of 
halo shapes as a function of redshift. Doing this it is convenient to parametrize a halo 
by its virial mass, M, and concentration parameter, c, instead of the parameters ps and 
Ts given in Equation ([2]). The concentration parameter is defined as c = -R/r_2, where 
r_2 is the distance where the profile falls as (r_2 = for the NEW profile) and R is 
the virial radius within which the halo has the mass M and a mean density of S{z)p{z), 
where p{z) again is the mean background matter density. The concentration parameter 
is dependent both on mass and redshift and we model this dependence according to 
|97] . Here the concentration parameter is also treated as a stochastic variable with a 
log- normal distribution for a fixed mass. Simplified, we can write the quantity A'^{z) 
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Figure 23. Spectral shapes of the EGRET measurements of the EGRB, with (51] 
(dots) and ,52j (squares) used for our conservative background, and 94J (dashed) for 
the the unresolved blazar model, and two examples of cosmological WIMPs, Equation 
(fTO]) (including the effect of substructures, see text), with masses of 70 and 200 GeV. 
The dotted WIMP spectra are calculated with the absorption from [96] . 



where dn/dM is the halo mass function and -P(c) the log-normal distribution. Clump- 
ing the DM into halos typically yields a boost of 10'* < A^(2; = 0) < 10^, depending 
on the choice of halo profile and the model of halo concentration parameters. Both 
choices contribute about a factor of ten each to the uncertainty in the normalization of 
the cosmological WIMP signal. The concentration parameter model, together with the 
optical depth, changes not only the normalization of the cosmological signal but also its 
shape. This can be compared to point-sources where the spectrum can be divided into a 
product of two factors: one determining the shape, resulting only from particle physics, 
and one determining the normalization, originating from astrophysics, i.e. the shape of 
the density profile. The largest contribution to A^(z) comes from small halos formed in 
an earlier, denser Universe. However, our understanding of halos at the low mass end 
is limited due to finite resolution of the N-body simulation. Therefore we have to use 
a cut-off mass, below which we do not trust our toy models for the halo concentration 
parameters. We put this cut-off at IQ^Mq. Lowering the cut-off might boost the signal 
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even further but will also introduce further uncertainties. Obviously, the observation 
of a significant signal could yield important information about the DM structure at all 
redshifts. 

As mentioned in section 14.31 there should exist smaller, bound halos that have sur- 
vived tidal stripping within larger halos [78]. Although not as massive as the primary 
halos, the substructure halos arise in higher density environments, which makes them 
denser than their parent halos. Also, the subhalos are tidally stripped from the outside 
inward, which further increases their concentration. The phenomenon of halos within 
halos seems to be a generic feature, since detailed simulations reveal substructures even 
within sub-halos |98] . We model the subhalo structure by imposing that a certain frac- 
tion of the parent halo mass is in substructures, and we associate higher concentration 
parameters with them. Also here we only consider halos down to masses of approxi- 
mately IO^Mq. 

Fast detector simulations are performed (using the tools described in section 12.21) for a 
generic model of WIMPs with masses ranging from 50 GeV to 250 GeV. A analysis 
was performed to obtain a sensitivity plot in < av > vs ttl^jmp- The WIMP signal is 
computed using the NFW profile and with and without the effect of substructures. For 
the substructures, we assume that they constitute 5% of the mass and have four times 
the concentration parameter of their parent halo. As background to our signal we use 
both an optimistic spectrum of unresolved blazars [M] and a conservative assumption 
that the background is the EGRB as measured by [51] . We include the residual charged 
particle background in the detector at a level of 10 % (30 %) of the respective extra- 
galactic diffuse gamma ray background above (below) 10 GeV. 

The result (see figure shows that GLAST is sensitive to total annihilation cross- 
sections of the order 10~^^ — 10"^^ cm^ s~^, depending on the presence or absence of 
sub-halos. It should be noted that if the dominant fraction of DM were indeed thermal 
WIMPs annihilating according to our simplified model, cosmology would, to first order, 
constrain the cross-section to be < av >~ 3 • 10~^^ cm^ s~^ in order to get the correct 
relic density, which would be within the reach of GLAST. 

It should be noted that the predictions for background spectra from astrophysical 
sources are very uncertain, especially at high energies. Our sensitivity estimate also 
optimistically assumes that we have a perfect analysis in our determination of the 
EGRB. For the estimate of the signal expectation (as mentioned earlier), the exact 
shape of the halo is here less important than in other analyses: when assuming a Moore 
profile the sensitivity improves by a factor of roughly ten. 
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Figure 24. Str exclusion curves for one year of GLAST simulated data. The lower 
edge of the shaded band corresponds to a background as predicted by the Blazar 
model presented by [94j, and the upper bound corresponds to the conservative case, 
where the background flux is given by the analysis of EGRET data [ST . NEW denotes 
computation of the WIMP signal using a NEW profile ^26J and NEW+subhalo includes 
the effect of having 5 % of the host halo mass in substructures with four times higher 
concentration parameter than the parent halo. Note that a Moore profile would lead 
to an improvement in sensitivity by about a factor of ten. 



5. GLAST sensitivity to specific particle physics models 
5.1. GLAST and mSUGRA 

In the preceding sections we have considered a generic model for DM that causes the 
7-ray signal. We only assumed a WIMP that gives a mono-energetic quark (or fermion) 
spectrum or line features. In this context the cross-section and mass of the WIMP are 
free parameters. However, there are several extensions to the Standard Model of Parti- 
cle Physics (SM) which predict a particle state which could constitute the WIMP. The 
most studied class of such models is super symmetry, in particular minimal supersym- 
metric extensions to the SM (MSSM) and its constrained version mSUGRA, in which 
the soft supersymmetry breaking terms derive from a high-energy supergravity theory 
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with a common scalar mass mo and a common gaugino mass mi/2 at the GUT scale 
(for details see |99l|100]). In figure [25] we show the set of mSUGRA and MSSM mod- 
els which pass all accelerator constraints and are consistent with WMAP data in the 
parameters in which we calculated the LAT sensitivity, i.e. the {< av >, my/ imp) plane. 
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Ncutialino Mass (CcVj 

Figure 25. MSSM and mSUGRA models in the < av >,mwiMP plane. The models 
included in these regions are consistent with accelerator constrains and WMAP data. 
The lines represent the 5 a sensitivity from the GC (upper) and the 5 a sensitivity 
from a Galactic halo analysis (lower) corresponding to the best and worst sensitivities 
estimated in this paper for a NFW profile. 

For comparison with other experiments (for example accelerators) mSUGRA sensitivi- 
ties are often presented in the plane (mo,mi/2), where the other free parameters (tan/3, 
Ao and sgn{fi)) are fixed. Here tan /3 denotes the ratio of the vacuum expectation values 
of the two neutral components of the SU(2) Higgs doublet, Aq is the proportionality 
factor between the supersymmetry breaking trilinear couplings and the Yukawa cou- 
plings, while fi is determined (up to a sign) by imposing the Electroweak Symmetry 
Breaking (EWSB) conditions at the weak scale. To give an example we consider the 
LAT sensitivity to a signal from the GC. The analysis is similar to the one presented in 
section except that for ease of computations we use an approximation of the effective 
area instead of using the detector simulations. In particular, we assume a total exposure 
of 3.7 ■ 10^° cm^ s, an angular resolution (at 10 GeV) of ~ 3 • 10~^ sr, and 5 years of data 
taking. The error in the sensitive area occupied in the model space introduced by this 
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approximation is estimated to be about 30 %. The result, at a 3cr confidence level, for 
tan/5 = 55, Aq = and sgn{fi) = +1 is shown in figure [261 It can be seen that for 
this particular class of models, which is neither particularly optimistic nor pessimistic, 
GLAST is sensitive to a significant fraction of the cosmologically relevant parameter 
spaces. 

tan(3=55, Aq=0, sgn(|i)=+1 
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Figure 26. GLAST reach (solid lines) in the mSUGRA parameter space for tan/3 = 
55, Aq = and sgn(fi) = +1. The region below the lines can be excluded at Scr 
level after 5 years of data taking. Solid lines represent the GLAST reach for NFW 
(lower line) and Moore (upper line) profiles, while dashed lines represent the neutralino 
isomass contours (expressed in GeV) . Dark shaded regions are parameter space regions 
excluded by either theoretical or experimental (accelerator) constraints. The upper left 
region is excluded because in that region the lightest stau (not the neutralino) is the 
LSP, the lower left region is excluded due the accelerator bounds on the Higgs boson 
masses, b S7, slepton and squark masses, etc., while the right lower region is excluded 
because for these parameters there would be no electroweak symmetry breaking. 



5.2. GLAST and UED 

Universal extra dimension (UED) theories [lOlj provide another natural candidate for 
DM. In UED theories all Standard Model fields are allowed to propagate in a higher- 
dimensional bulk. After compactification of the internal space, these additional degrees 
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of freedom appear as towers of new, heavy states in the effective four- dimensional theory. 
The hghtest of these Kaluza-Klein particles (LKP) is usually stable due to KK parity, 
an internal symmetry analogous to R parity in the supersymmetry. Since for UED 
the chiral suppression of annihilation to light fermions is not present, significant yields 
(~ 20%) of electrons can be expected. We therefore study the KKDM case in the context 
of the LAT's capability to detect electrons (positrons). Based on the simulations that 
Baltz and Hooper |102j made using a NFW DM distribution with boost factor of 5 
and piocai=0.4: GeV /cm^ we can estimate the magnitude of the electron flux (after 
propagation) generated by the LKP annihilation (with the mass of rriLKp) as 



The flux produced by LKP annihilation would have a sharp cutoff in the electron spec- 
trum at the energy corresponding to the LKP mass. The observational goal would be to 
see a signal with a sharp cutoff riding atop the "conventional" electron spectrum. With 
the large LAT effective area, the number of electrons detected could be sufficient to 
provide a statistically detection of such a spectral feature. Figure [27] shows the compar- 
ison of that "conventional" electron flux [103j , [104\ with the flux from LKP annihilation 
(peak value). 
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Figure 27. Expected electron flux from LKP annihilation, along with the observed 



electron flux (see [103] and references therein). 



Using a simple sequential cut analysis using topological variables provided by the tracker 
and the calorimeter, it is possible to detect cosmic ray electrons efficiently [104] . In the 
energy range from 20 GeV to ~ 1 TeV the LAT effective geometric factor for electrons is 
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0.2 — 2ni?sr, and the energy resolution is 5-20%, depending on the energy. The residual 
hadron contamination in this analysis is ~ 3% of the remaining electron flux. The 
expected number of electrons is ~ 10'' electrons per year in 20 GeV - 1 TeV energy 
range. The observation time needed to detect an LKP feature within the assumed 
model with Sa signiflcance is shown in flgure [281 An LKP with a mass of ~ 600 GeV 
is probably the heaviest one which can be observed within the assumed constraints and 
the GLAST mission duration. Taking into account that for thermal freeze-out [102] 
the preferred LKP mass is in the range 600-700 GeV, the window for detecting a LKP 
mass in the LAT search is not large. However, similar analyses of the LAT capability to 
detect DM annihilation can be applied to any model in which electrons are produced. 




400 500 600 

LKP mass, GeV 

Figure 28. LKP detection in the LAT electron spectrum. The time needed to detect 
an LKP feature with 5cr significance. 

As an illustration, we simulate the electron spectrum to be detected by the LAT in 1 
year of observation 10^ electrons). For this illustration we consider a signal from the 
nearest DM clump at a distance of 100 pc, taking into account diffusive propagation of 
electrons [105] and masses of the LKP of 300 and 600 GeV. The signal can be identifled 
with high statistical signiflcance on top of the "conventional" electron flux [103] . p^^ 
(as can be seen in flgure flgure [29]) . This is a very favorable model, but recent results of 
the balloon-borne instruments ATIC [106] and PPB-BETS [107] indicate the presence 
of such a feature in the electron spectrum at 300-500 GeV. 
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Figure 29. Simulated detection of LKPs with masses of 300 GeV and 600 GeV in 
the LAT electron spectrum to be collected in 5 years of operation. Filled circles — 
"conventional" electron flux; open circles — the same but with added signal from 300 
GeV LKP, and open squares — the same with added signal from 600 GeV LKP. For 
this signal, only the nearest clump at a distance of 100 pc is considered. The breaks 
are also shown by dashed lines to guide the eye. 



6. Discussion & Conclusions 

Using the current state-of-the-art Monte Carlo and event reconstruction software de- 
veloped within the LAT collaboration, we present preliminary sensitivity calculations 
for several astrophysical searches of DM annihilation. In particular, we present sensi- 
tivities for detecting DM annihilation from the GC, Galactic and extragalactic diffuse 
emission. Galactic known and unknown satellites, point sources and a dedicated search 
for the 27 line signal. We exemplify the possibility to constrain specific particle physics 
scenarios (especially within in mSUGRA and UED frameworks) on the search for DM 
annihilation 7-ray emission from the GC and by using the LAT not as a 7-ray detector, 
but employing its very good electron/positron detection capabilities at very high energy. 

The sensitivities presented here are based on analyses which are idealized in the sense 
that systematic uncertainties in the instrument performance estimates are neglected and 
that uncertainties in the background estimate are quantitatively taken into account only 
in a preliminary manner. For a given particle physics model, the dominant uncertainty 
in the sensitivities presented here is due to the lack of knowledge on the DM density 
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distribution. Depending on which DM halo profile and/or substructure is assumed sen- 
sitivities can easily vary by one or two orders of magnitude. For most examples we 
consider "standard" DM halo structure, i.e. we are conservative in the sense that we 
do not include density enhancements that might boost the expected annihilation signals. 

Using simplifying assumptions, the uncertainty in the background prediction (most cru- 
cial in the search for DM signal from the halo) leads to a decrease in sensitivity by 
between 25 % and 45 % depending on the mass of the WIMP. For the GC analysis, 
the dominant background will be from sources in the vicinity of the GC, removal of 
which will lead to systematic uncertainties. Assessing those without the GLAST data 
does not make much sense, given uncertainties in extrapolating into the GLAST energy 
region. For the diffuse extragalactic and high-latitude source searches the charged par- 
ticle background and uncertainties therein are potentially important. Prior to launch, 
the levels of the charged particle background and its uncertainty are very difficult to 
assess. Estimates have to rely on poorly constrained Monte Carlo simulations. The 
charged particle background included at the level of roughly 10 % of the extragalactic 
background, which is in compliance with specifications, leads to about a 10 % decrease 
in sensitivity for a signal in the EGRB and in the number of detectable satellites. Un- 
certainties in the charged particle background are negligible for the sensitivity for both 
signals (assuming they are at the level of ~ 20 %). 

For the Galactic and extragalactic diffuse background the range of possible backgrounds 
is illustrated by assuming several representative models, which are discussed in detail. 
The sensitivities are preliminary in the sense that estimates of the instrument perfor- 
mance, analysis methods and estimates of the expected backgrounds are being contin- 
uously improved. In context of specific particle physics scenarios, also the estimates 
of the signals are continuously updated: For example, calculations of the 7-ray flux for 
Supersymmetric Dark Matter annihilation incorporating QED corrections, indicate that 
for part of the parameter space 7-ray yields might be boosted by three or four orders of 
magnitude and and lead to distinct spectral signatures [108j . 

We conclude that the LAT has good potential to discover DM annihilation for a signifi- 
cant fraction of interesting parameter space, i.e. for values of annihilation cross-section 
of between < av >^ 10^^^ cm^s~^ and < av >^ 10~^^ cm^s~^ depending on WIMP 
masses in the range between 40 and 500 GeV. For less conservative assumptions on 
the Dark Matter density (for example additional substructure or a Moore profile) the 
sensitivity improves by one to two order of magnitudes. Our conclusions are consis- 
tent with previous work that employed cruder representations of the GLAST response 
and/or less thorough treatments of the backgrounds (see e.g [6], [IE], [69] and references 
therein). If indeed a significant DM signal is present, GLAST will be able to image the 
DM structure in our Galaxy. 



46 



Acknowledgments 

We thank the members of the LAT collaboration for providing many interesting dis- 
cussions and elements of the software used in this analysis. Sergio Colafrancesco, Steve 
Ritz and Julie McEnery are thanked for careful reading of the manuscript. Numerous 
discussions with Jeff Scargle are gratefully acknowledged. This work was supported by 
the U.S. Department of Energy contract number DE-AC02-76SF00515, the U.S. De- 
partment of Energy contract number, DE-FG02-91ER40690, NASA, Vetenskapsradet 
and the Swedish Space Board. IM acknowledges partial support by the NASA APRA 
program. The LAT is being developed by an international collaboration with primary 
hardware and software responsibilities at (in alphabetical order) The Agenzia Spaziale 
Italiana, Centre National de la Recherche Scientifique / Institut National de Physique 
Nucleaire et de Physique des Particules, Commissariat a I'Energie Atomique, God- 
dard Space Flight Center, Hiroshima University, Istituto Nazionale di Astrofisica, Isti- 
tuto Nazionale di Fisica Nucleare, Naval Research Laboratory, Ohio State University, 
Kalmar University, Royal Institute of Technology - Stockholm, Stanford Linear Acceler- 
ator Center, Stanford University, University of California at Santa Cruz, and University 
of Washington. The LAT project is managed by the Stanford Linear Accelerator Cen- 
ter, which is also the location of the Instrument Science Operations Center, and the 
LAT Principal Investigator is Peter Michelson at Stanford University. Other institutions 
that have made significant contributions to the instrument development include the In- 
stitute of Space and Astronautical Science, Stockholm University, University of Tokyo, 
and Tokyo Institute of Science and Technology. 

References 

[1] W. B. Atwood [GLAST Collaboration], Nucl. Instrum. Metli. A 342, 302 (1994). 
[2] P. F. Michelson et ai, GLAST Proposal to NASA, Response to AO 99-OSS-03, Stanford University, 
(1999). 

[3] P. F. Michelson [GLAST-LAT Collaboration], AIP Conf. Proc. 921 (2007) 8. 
[4] C. Meegan et al, AIP Conf. Proc. 921 (2007) 13. 
[5] |http : //www-glast . staf ord . edu| 

[6] G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) 

PXiv:hep-ph/9506380|; L. Bergstrom, Rept. Prog. Phys. 63, 793 (2000) 
Xiv:hep-ph/0002126|; G. Bertone, D. Hooper and J. Silk, Phys. Rept. 405, 279 (2005) 
arXiv:hep-ph/0404175|. 
[7] D. Hooper and S. Profumo, Phys. Rept. 453, 29 (2007) |a rXiv:hep-ph/0701197| . 
[8] http://cossc.gsfc.nasa.gov/docs/cgro/cgro/egret.html 

[9] E. A. Baltz, M. Battaglia, M. E. Peskin and T. Wizansky, Phys. Rev. D 74, 103521 (2006) 
arXiv:hep-ph/0602 187] . 

[10] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 377 (2007) 
[arXiv: astro-ph/0 603449| . 

[11] S. Perlmutter et al. [Supernova Cosmology Project Collaboration], Astrophys. J. 517, 565 (1999) 
arXiv:astro-ph/9812133|; A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 
116, 1009 (1998) arXiv:astro-ph/9805201] . 



47 



M. Tegmark et ai, Phys. Rev. D 74, 123507 (2006) |arXiv:astro-ph/0608632 



A. G. Sanchez et ai, Mon. Not. Roy. Astron. See. 366, 189 (2006) arXiv:a stro-ph/0507583| . 
M. Bradac et ai, Astrophys. J. 652, 937 (2006) arXiv:astro-ph/0608408]. 



S. Hannestad, JCAP 0305, 004 (2003) |arXiv:astro-ph/0303076 



L. Bergstrom, P. UUio and J. H. Buckley, Astropart. Phys. 9, 137 (1998) 'arXiv:astro-ph/9712318]. 
M. Gustafsson, E. Lundstrom, L. Bergstrom and J. Edsjo, Phys. Rev. Lett. 99 041301 (2007) 

arXiv:astro-ph/0 703512| . 
A. A. Moiseev et al , Astropart. Phys. 27 (2007) 339 jarXiv:astro-ph/070258T] . 
http:/ /glast. gsfc.nasa.gov/ssc/data/ 



http://www-glast.slac.stanford.edu/software/IS/glastJat_perforinance.htin 



B. C. Allanach, S. Kraml and W. Porod, JHEP 0303 (2003) 016 [arXiv :hep-ph/0302102|. 



G. Belanger, S. Kraml and A. Pukhov, Phys. Rev. D 72, 015003 (2005) |arXiv:hep-ph/0502079l. 
P. Gondolo, J. Edsjo, P. UUio, L. Bergstrom, M. Schelke and E. A. Baltz, JCAP 0407 (2004) 008 
arXiv:astro-ph/0406204| . 

A. Cesarini, F. Fucito, A. Lionetto, A. MorseUi and P. Ulho, Astropart. Phys. 21 (2004) 267 
Pythia program package, see T. Sjostrand, Comp. Phys. Comm. 82 (1994) 74. 

J. F. Navarro, C. S. Frenk and S. D. White, Astrophys. J. 462 (1996) 563 

B. Moore, T. Quinn, F. Govcrnato, J. Stadel and G. Lake, Mon. Not. Roy. Astron. Soc. 310 (1999) 
1147 

W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 
A. W. Strong, Astrophys. Space Sci. 309 (2007) 35 |arXiv:astro-ph/0609359). 
A. W. Strong, L V. Moskalenko and O. Reimer, Astrophys. J. 537 (2000) 763 
A. W. Strong, L V. Moskalenko and O. Reimer, Astrophys. J. 613, 962 (2004) 
'arXiv:astro-ph/0406254]. 



Galprop,2008, http://galprop.stanford.edul 

A. W. Strong, I. V. Moskalenko and V. S. Ptuskin, Ann. Rev. Nucl. Part. Sci. 57, 285 (2007) 
arXiv: astro-ph/070 1 5 1 7| . 

L V. Moskalenko, A. W. Strong, and O. Reimer Diffuse Gamma Rays: Galactic and Extragalactic 
Diffuse Emission ed K S Cheng and G E Romero {Dordrecht: Kluwer) Astrophysics and Space 
Science Library vol 304 (279-310) p 2004 
S. D. Hunter et al, Astrophys. J. 481, 205 (1997). 

L V. Moskalenko, S. W. Digel, T. A. Porter, O. Reimer and A. W. Strong, Nucl. Phys. Proc. 

Suppl. 173 (2007) 44 [arXiv:astro-ph/0609768| . 
F. W. Stecker, S. D. Hunter and D. A. Kniffen, arXiv:0705.4311 [astro-ph]. 

B. M. Baughman, W. B. Atwood, R. P. Johnson, T. A. Porter and M. Ziegler, larXiv:0706."0503l 
[astro-ph]. , Contributed to 30th International Cosmic Ray Conference (ICRC 2007), Merida, 
Yucatan, Mexico, 3-11 Jul 2007. 

F. Yiou et al 1997, J. Geophys. Res.-Oceans, 102, 26783 
K. Knic et al 2004, Phys. Rev. Lett., 93, 171103 

G. Case and D. Bhattacharya 1996, A&A Suppl., 120C, 437 

L V. Moskalenko and A. W. Strong, AIP Conf. Proc. 801 (2005) 57 arXiv:astro-ph/0509414]. 
A. W. Strong and L V. Moskalenko, Astrophys. J. 509, 212 (1998) arXiv:astro-ph/9807150] . 
E. S. Seo and V. S. Ptuskin 1994, ApJ, 431, 705 

S. Orito et al. [BESS Collaboration], Phys. Rev. Lett. 84, 1078 (2000) [arXiv:as tro-ph/990642^. 
L V. Moskalenko, A. W. Strong, J. F. Ormes and M. S. Potgieter, Astrophys. J. 565, 280 (2002) 

arXiv:astro-ph/0106567|. 
T. Hams et al. [BESS Collaboration], Proc. 30th Int. Cosmic Ray Conf. (Merida) 
I. V. Moskalenko, A. W. Strong, S. G. Mashnik and J. F. Ormes, Astrophys. J. 586, 1050 (2003) 
arXiv: astro-ph/02 104801 . 

[49] I. V. Moskalenko, A. W. Strong, S. G. Mashnik and J. F. Ormes, Proc. 28th Int. Cosmic Ray 
Conf. (Tsukuba) 4 1921 2003 |arXiv:astro-ph/0366368[ 



48 



D. J. Thompson and C. E. Fichtel 1982, A&A, 109, 352 

P. Sreekumar et al [EGRET Collaboration], Astrophys. J. 494, 523 (1998) 

arXiv:astro-ph/9709257] . 
A. W. Strong, I. V. Moskalenko and O. Reimer, Astrophys. J. 613, 956 (2004) 

arXiv: astro-ph/040544l| . 
I. V. Moskalenko and A. W. Strong, Astrophys. J. 528, 357 (2000) 'arXiv:astro-ph/9811284l. 
I. V. Moskalenko, T. A. Porter and S. W. Digel, Astrophys. J. 652, L65 (2006) [Erratum-ibid. 
664, L143 (2007)] |arXiv:astro-ph/060752l] . 

E. Orlando and A. Strong, Astrophys. Space Sci. 309, 359 (2007) | arXiv:astro-ph/0607563 . 

E. Orlando, D. Petry and A. Strong, AIP Conf. Proc. 921, 502 (2007) [ arXiv:0 704.0462 [astro-ph]]. 
I. V. Moskalenko et al. larXiv:0712.2015 [astro-ph]. 

F. W. Stacker and M. H. Salamon, Astrophys. J. 464, 600 (1996) arXiv:astro-ph/9601126 . 
D. Elsaesser and K. Mannheim, Phys. Rev. Lett. 94, 171302 (2005) ar Xiv:astro-ph/0405235i . 
W. de Boer, New Astron. Rev. 49, 213 (2005) arXiv:hep-ph/0408166 . 

F. Aharonian et al (HESS CoUaboration) 2004, A&A, 425, L13 
J. Albert et al (MAGIC Collaboration) 2006, ApJ, 638, LlOl 

K. Kosack, et al (VERITAS CoUaboration) 2004, ApJ, 608, L97 

T. Tsuchiya, et al (CANGAROO Collaboration) 2004, ApJ, 606, L115 

S. Profumo, Phys. Rev. D 72 (2005) 103521 ar^v^^stro-ph/ 0508628]. 

L. Bergstrom, T. Bringmann, M. Eriksson and M. Gustafsson, Phys. Rev. Lett. 94, 131301 (2005) 
arXiv:astro-ph/04 10359| . 

C. van Eldik, O. Bolz, I. Braun, G. Hermann, J. Hinton and W. Hofmann (H.E.S.S. Collaboration), 
arXiv:0709.3729 [astro-ph]. Contributed to 30th International Cosmic Ray Conference (ICRC 
2007), Merida, Yucatan, Mexico, 3-11 Jul 2007. 

D. Hooper and B. L. Dingus, Phys. Rev. D 70 (2004) 113007 |arXiv:astro-ph/021061'^ . 
S. Dodelson, D. Hooper and P. D. Serpico. [arXiv:0711.4621l [astro-ph]. 

H. Mayer-Hasselwandcr et al , Astron. Astrophys. 335 (1998) 161. 
P. Gondolo et al , JCAP, 0407 (2004) 008 

D. Hooper, G. Zaharijas, D. P. Finkbeiner and G. Dobler. lar"Xiv:0709.3lTl [astro-ph]. 

D. P. Finkbeiner, Astrophys. J. 614 (2004) 186 [arXiv:astro-ph/031154';^ . 
M. Regis and P. UUio, arXiv:0802.0234 [hep-ph]. 

G. R. Blumenthal, S. M. Faber, J. R. Primack and M. J. Rees, Nature 311, 517 (1984). 
P. J. E. Peebles, Astrophys. J. 284, 439 (1984). 

S. Profumo, K. Sigurdson and M. Kamionkowski, Phys. Rev. Lett. 97 (2006) 031301 
arXiv:astro-ph/0603373 | . 

E. Hayashi, J. F. Navarro, J. E. Taylor, J. Stadel and T. Quinn, Astrophys. J. 584 (2003) 541 
arXiv:asto-ph/0203004| . 

J. E. Taylor and A. Babul, Mon. Not. Roy. Astron. Soc. 364, 535 (2005) arXiv:astro-ph/0410049]. 

J. E. Taylor and A. Babul, Mon. Not. Roy. Astron. Soc. 364, 515 (2005) 

|arXiv:astro-ph/0410048| . 
V. R. Eke, J. F. Navarro and M. Steinmetz, Astrophys. J. 554, 114 (2001) 

iarXiv:astro-ph/0012 337| . 
J. Diemand, B. Moore and J. Stadel, Nature 433, 389 (2005) [arXiv:astro-ph/050158^ . 
G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267, 195 (1996) 

arXiv:hep-ph/9506380l . 
L. Fieri, G. Bertone and E. Branchini. larXiv:0706. 21011 [astro-ph]. 
J. Diemand, M. Kuhlen and P. Madau, Astrophys. J. 657 (2007) 262 
J. E. Taylor and J. Silk, Mon. Not. Roy. Astron. Soc. 339 (2003) 505 
E. Berlin and S. Arnouts, S., 1996, 



arXiv:astro-ph/0611370 



arXiv:astro-ph/0207299 



A&AS, 117, 393. See also http://terapix.iap. fr/rubrique.php?id_rubrique=91/ 



[87] E. A. Baltz, J. E. Taylor and L. L. Wai, arXiv:astro-ph/0610731 



49 



[88] N. W. Evans, F. Ferrer and S. Sarkar, Phys. Rev. D 69, 123501 (2004) [arXiv:astro-ph/0311 145] . 
[89] R. C. Hartman et al. [EGRET Collaboration], Astrophys. J. Suppl. 123 (1999) 79. 
[90] J. T. Kleyna, M. I. Wilkinson, N. W. Evans and G. Gilmore Astrophys. J. 563, L115 (2001) 
|arXiv:astro-ph/0111329 |. 

[91] R. A. Ibata, R. F. G. Wyse, G. Gilmore, M. J. Irwin, N. B. Suntzeff, Astron. J. 113, 634 (1997) 
arXiv:astro-pli/9612025|. 

[92] G. Bertone, A. R. Zentner and J. Silk, Phys. Rev. D 72, 103517 (2005) [arXiv:astro-ph/ 0509565] . 
[93] F. Stoehr, S. D. M. White, V. Springel, G. Tormen and N. Yoshida, Men. Not. Roy. Astron. Soc. 

345 (2003) 1313 [arXiv:astro-ph/0307026|. 
[94] P. UUio, L. Bergstrom, J. Edsjo and C. G. Lacey, Phys. Rev. D 66, 123502 (2002) 

[arXiv:astro-ph/0207125|. 

[95] J. R. Primack, R. S. Somerville, J. S. Bullock and J. E. G. Devriendt, AIP Conf. Proc. 558, 463 
(2001) arXiv: astro-ph/00 1 1475| . 



F. W. Sleeker, AIP Conf. Proc. 921, 237 (2007) arXiv:astro-ph/0703505 



[97] J. S. Bullock et al , Mon. Not. Roy Astron. Soc. 321 (2001) 559 arXiv:astro-ph /9908159 



[98] J. Diemand, M. Kuhlen and P. Madau, Astrophys. J. 667, 859 (2007) [arXi v:astro-ph/070333^ . 
[99] L. J. HaU, J. D. Lykken and S. Weinberg, Phys. Rev. D 27, 2359 (1983). 
[100] N. Ohta, Prog. Thcor. Phys. 70, 542 (1983). 

[101] T. Appelquist, H. C. Cheng and B. A. Dobrescu, Phys. Rev. D 64 (2001) 035002 

arXiv:hep-ph/0012100l . 
[102] E. A. Baltz and D. Hooper, JCAP 0507, 001 (2005) [arXiv:hep-ph/0411053| . 
[103] S. Torii et al, Astrophys. J. 559 (2001) 973. 

[104] A. A. Moiseev, J. F. Ormes and I. V. Moskalenko, Proc. 30th Int. Cosmic Ray Conf. (Merida) 

arXiv:0706.088 2 [astro-ph]. 
[105] A. A. Moiseev et al , Frascati Phys.Ser.45:l-328,2007. Edited by A. Lionetto, A. Morselh ( ISBN 

978-88-86409-54-0) 
[106] J. Chang et al , 2005, Proc. 30th Int. Cosmic Ray Conf. (Pune) 
[107] K. Yoshida et al , 2005, Proc. 30th Int. Cosmic Ray Conf. (Pune) 
[108] T. Bringmann, L. Bergstrom and J. Edsio. [arXlv:0710.3169l [hep-ph]. 

[109] H. C. Cheng, J. L. Feng and K. T. Matchev, Phys. Rev. Lett. 89 (2002) 211301 
.arXiv:hep-ph/0207125j . 



